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o: 

. The explicit solution to the spectral problem of quantum graphs is used to obtain 

(N 

(3jT)! the exact distributions of several spectral statistics, such as the oscillations of the 

quantum momentum eigenvalues around the average, 5k n = k n — k n , and the nearest 
Q\ ■ neighbor separations, s n = k n — fc n _i. 
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i I. INTRODUCTION 

o 

Understanding statistical properties of the quantum spectra produced by the classically 



nonintegrable systems is one of the most fundamental problems of quantum chaos theory. 



According to the Random Matrix Theory (RMT), the statistics of the spectral fluctuations is 
universal, and reflects only the general symmetry properties of system The universality 



of the RMT approach and its ability to produce specific predictions are particularly valuable 
since the spectra of the systems to which it is applied are typically out of reach for direct 
analytical studies. On the other hand, understanding how the RMT distributions themselves 
emerge from the phase space structures used by semiclassical spectral theories, such as 
the periodic orbit theory, is a matter of intense research. There has been a number of 
publications Q] dedicated to computing spectral statistics that are accessible via periodic 



orbit expansion for t 



le density of states. This analysis was particularly complete for the 



quantum graphs , which are simple and convenient models of quantum chaos. 

As a reminder, quantum graphs consist of a quantum particle moving on a quasi one- 
dimensional network with B bonds and V vertexes. In the limit h = 0, these systems produce 
a nonintegrable (mixing) classical counterpart - a classical particle moving on the same net- 
work, scattering randomly on its vertexes. Despite the apparent simplicity such system, its 
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classical behavior exhibits many familiar features of multidimensional deterministic chaotic 
systems jj], which is clearly manifested in quantum regime. Numerical investigations of the 
quantum graph spectra p| show that the spectral distributions of sufficiently complicated 
networks are closely approximated by the RMT Wignerian distributions. On the other 
hand, periodic orbit theory for quantum graphs yields exact harmonic expansions for the 
density of states, spectral staircase, quantum and classical zeta functions, etc. Moreover, as 
shown recently in ^], the periodic orbit theory for the quantum graphs can be "localized" 
- it is possible to obtain the harmonic series expansion representation for the individual 
eigenvalues of the energy or of the momentum, k n = k(n), as a global function of the index 
n. Hence, despite classical nonintegrability the spectral problem for quantum networks is 
exactly solvable within the periodic orbit theory approach. This fact provides an interesting 
opportunity to obtain analytically several new spectral statistics in terms of the periodic 
orbit theory, which is the main subject of this Letter. 

II. SPECTRAL DISTRIBUTIONS FOR REGULAR QUANTUM GRAPHS 

In the simplest case of the so called regular graphs 0, [?J the individual momentum 
eigenvalues can be expanded into a periodic orbit series, 

k n = j-n - ^- — sin ( ^f) sin (u p n) . (1) 



Lq Lq uj p \ 2 



Here L is the total action length of the graph, A p is the weight factor o: 



which is a function of the scattering coefficients at the graph vertexes . The frequency 
uj p = 7r {rripLi + ... + nip L B ) /Ln, is defined by the numbers of times, m p , the orbit p passes 
over the bond i. As shown in p, |6(, the series in is convergent and bounded by vr/Lo, 
so the values k n are locked within a sequence of periodic cells. The first term in (^J) gives 
the average (Weyl) behavior of the eigenvalue k n and the subsequent periodic orbit sum 
describes the fluctuations 5n^ = L (k n — k n ) /ir around the average. 

The transition from the exact formula (0) to statistical description of Sn^s can be made 
based on the well known fact from the analytical number theory |9( that the sequence of the 
remainders x n = [cm| modl , n — 1, 2, for any irrational number - a is uniformly distributed 
over the interval x G [0,1]. Assuming the generic case in which every f2j — L j / Lq is an 
irrational number, it is clear that parsing through the spectral sequence k±, k 2 , k n ,..., will 



a periodic orbit p, 
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FIG. 1: The histogram of the spectral fluctuations for the simplest regular two bond star graph 
obtained from the exact expansion Q (shown in red), compared to the distribution obtained from 
(j3J) using the orbits that include up to 6 scattering events (shown in blue). 

generate a sequence of random phases in ([!}, denned via the combinations 

Xi = [7ran] mod27r , (2) 

that are uniformly distributed in the interval [0, 2ir]. Hence, the deviations of the eigenvalues 
from the average, 5n\ are statistically described by a series of random inputs 

5 ( 0) = _2^A, si /^ sin 

71 * — ' UO v V L I 

V ' 



where x = (xi, ...,xb), which has the same structure as the periodic orbit sum (JTJ) (see 
Fig. 1). The distribution of <5's can then be obtained as 

pf (5W) = J 2 J § ^(0) + Cf sin (m p x) j dx (4) 

where = l^sin^ and dx = Yli^- Using the exponential representation of the 5 - 
functional, one finds that 

pV>(8®) = Jdke m(0) F^(k), (5) 



where 

»2tt 



F^\k) = I e^ c ^ s ' m ^T\dx (6) 
Jo 
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FIG. 2: The histogram of the nearest neighbor separations for the regular two bond star graph 
obtained via direct numerical solution of the spectral equation (shown in red), compared to the 
distribution obtained from (|SJ) (shown in blue). To emphasize the connection with Fig. 1, the graph 
is scaled and shifted along the abscissa by tt/Lq. 

is the characteristic function of the distribution. As in the case of the series expansion for 
k n , a finite order (Ith) approximation to the exact result in (jHJ) is obtained by considering 
only the orbits that involve a particular number (/) of vertex scatterings jf],!?!]. 

Similarly, for the spacings between two eigenvalues, = k n+m — k n , the periodic orbit 
expansion, 

4% = f rn - J2 Dfl cos (u, p (» - |) ) , (7) 
v 

where -Dp°m = -^^j-sin (^y) sin (^^), yields the distribution Pf2 of the type where 
now 




Note, that despite the complete nonintegrability of the underlying classical system |4| , the 
statistics of the spectral fluctuations of certain (e.g. some regular) quantum graphs is not 
described by the RMT and so the shapes of their probability distributions deviate noticeably 
from the familiar Wignerian profiles (Fig. 2). 

All the above distributions are exact and depend explicitly on the specific graph param- 
eters - the vertex scattering coefficients and the periodic orbit set. 
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III. SPECTRAL HIERARCHY METHOD 

As shown in |j. a generic quantum graph systems is irregular. The idea of producing the 
statistical distribution profiles for irregular quantum graphs is based on three key properties 
of the spectral determinant A (A;) [5], briefly outlined below. First, its roots as well as the 
roots of all of its derivatives are real. Secondly, there is exactly one root of A^\k) between 
every two neighboring roots of A^ +1 ^(/c) |8|. This implies that the zeroes of A(k) can be 
bootstrapped by the zeroes of A'(fc), while the latter sequence can be bootstrapped by 
the roots of A"(k) and so on (Fig. 3). Hence, there exists a hierarchy of almost periodic 
sequences, k { J ] , j = 0, 1, A^ (k%A = 0, such that fcjli < k { J~ 1] < k ( J ] and 

W = + 0) 

Lastly, the higher is the order j of the derivative, the smaller are the fluctuations of S n ^- 
In fact, the roots of a sufficiently high order r of the derivative, A^(fc), can be locked 
between a periodic bounding sequence of points, just as the spectra of the regular graphs 
sj]. Hence after r steps there is no need to consider the roots of A ( - r+1 \k) to separate kn^ 
from one another. Instead, one can use the periodic sequence ©with^ r+1) = ^(n+|), 
which explicitly carries the index n. The regular graphs discussed in the previous section 
correspond to the case r = 0. These three properties of the auxiliary sequences kn allow 
to obtain the whole hierarchy explicitly. Starting from the periodic separators, one can find 

* (r) 

the roots kn of the rth derivative of spectral determinant, then use them as separators on 
the next level of the hierarchy to find k n r , and so on. After r + 1 steps, 

h$- x) = / p ij - 1] (k)kdk, (10) 

n 

j — r + 1, 1, the spectrum k n is produced p. The integrations (fTUj) are made explicit 
by using the series expansions for the density p^\k) of the separators on each level of the 
hierarchy, 

p U) {k) = h + Re Y^S^A^e^ )k . (11) 

p 

The weight coefficients Ap and expansion frequencies in the expansions (jllj) are different 
for the different levels j of the hierarchy For j = the series ()11|) is just the usual Gutzwiller's 
formula. 
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FIG. 3: The bootstrapping of the spectral staircase (shown in red) by the N^\k) staircase (shown 
in green), bootstrapped in turn by the N^ 2 \k) staircase (shown in blue). Small vertical shift of 
N^{k) and N&\k) is given for illustration purposes. Note that the N^(k) does not bootstrap 
the N^°\k). Data obtained for the "fully connected quadrangle" graph 

Using © and (fTTf in one can define each element in the jth separating sequence in 
terms of the elements of the previous one. The result is a set of equations that completely 
describe the propagation of the fluctuations across the hierarchy, 

5t 1] = ft 1] - E C t 1] ™ + , (12) 

p 

where now the "zero frequency" term, /} J ~ 1) = [6$ - 8^l_ x - (Sn^) 2 + (#n2i) 2 J / 2 , the ex- 
pansion coefficients, Cp" 1 ^ = ^^f—jy sin ^ — — <^f-i + ^j-> anc ^ the phases, 

Vt 1] = + - l) /2 (13) 

are functions of the fluctuations 5„ and on the previous level of the hierarchy. In the 
particular case when r = 0, 5„ = 1/2, n = 1, (fT2"|) coincides with the oscillating part of 
©• 

The properties of the fluctuation sequences outlined in this Section are inherited by 
other spectral characteristics, such as the nearest neighbor separations (Fig. 4), etc. The 
corresponding systems of equations that describe the propagation of these characteristics 
across the hierarchy are similar to ()12j) and can be easily derived from relating a particular 
spectral sequence to the one that precedes it in the hierarchy via (fTUj) . 
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FIG. 4: The development of the full scale spectral fluctuations for the nearest neighbor separations 
(s(°), shown in red) for the fully connected quadrangle graph studied also in jjj for r = 6. One can 
follow the appearance of the Wignerian type distribution at the physical level from the distributions 
P^(s). 

IV. FLUCTUATION STATISTICS 

The expansions ()12j) allow to produce the probability distributions for the spectral fluc- 
tuations across the hierarchy, using the same approach as in Section II. By construction, the 
index n is the same across the hierarchy, so the combinations (J2J) in the arguments of each of 
the trigonometric expansion terms in ()12j) produce uniformly distributed random variables 
X{. This yields a discretized Ito type equation 

8^ = ft 1] ~ E CM sin (m^x + <pW) , (14) 
p 

where 8[ and Slf are also functions of x, produced by the corresponding expansions for 
and 5„2i- For a given j, these functions introduce the fluctuations from the lower levels of 
the hierarchy into ()14j) . 

It is possible to use the systems of expansions (|T^|) directly to produce the probability 
distributions for each of the S^'s. However, it is more illustrative to use instead a simple 
physical approximation, in which the jth level fluctuations 5n and are treated as 
independent random variables 8^ and distributed according to Pg. This allows to 
follow the accumulation of the fluctuations of the final distributions for 8^°' starting from 
the fluctuations at the regular level, S^ r \ At the regular level, j = r, the fluctuations 8^ 
are described by an expansion structurally identical to Q so the corresponding distribution 
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P 5 is self-contained and has the same functional form as (0). Knowing Pg allows to 
obtain pj?~ x \ and so on, up to the distribution P$°\ which applies to the fluctuations of 
the physical momentum eigenvalues. With the independent 8± and 8^ the transition from 
P^ to P^ is direct. Proceeding as in (JHJ) -© yields the result 

Pt 1 ^) = J dke^ (Ft^ik))^ , (15) 

where Fg'~ 1 \k,8^) is obtained as © for the series (JHJ), and (*)qo) represents averaging 
over the "disorder" produced by the separators 8^ and 8% of jth level using the weight 
function 

n c--D = Hti-vyp^pV) ^ P u) ^ (16) 

Similar considerations produce the distributions for the nearest neighbor separations, the 
form factor, etc. 



V. DISCUSSION 



The explicit periodic orbit representations of the individual momentum eigenvalues k n of 
quantum graphs [5] provide a regular analytical method of describing the statistical proper- 
ties of their spectra within the standard periodic orbit theory framework. 

As shown in j^j, solving the spectral problem for a generic quantum graph goes beyond 
the conventional Gutzwiller's theory approach. The additional information supplied by the 
auxiliary separating sequences kn allows to produce the exact physical spectrum k n = kn\ 
to reveal the hierarchical structure of the spectral fluctuations and then to unfold the exact 
probability distributions for spectral statistics associated with each sequence kii ^ ■ 

Among other things, the expansions (J3J) also provide a physical understanding about the 
origins of the universal behavior of spectral characteristics based on the periodic orbit the- 
ory. It is known that the terms of lacunary trigonometric series behave as weakly dependent 
random variables (see e.g. U] and the references therein). This fact allows to establish 



a number of universal features for the asymptotic distributions of their sums, such as conver- 
gence to Gaussian distribution with a specific variance (also conjectured to be the universal 
)robability distribution profile for the distribution of generic quantum chaotic systems 



12|). In addition, the build up of the distributions PW' along the levels of the hierarchy can 
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lead to the appearance of other universal, (albeit more complex, e.g. Wignerian, see Fig. 4 
and l|) profiles. 

The advantage of this result is twofold. On the one hand, since the spectral properties 
of sufficiently complex quantum graphs are well described by the RMT j^j], the proposed 
approach provides a mechanism to follow the build up of the universal distributions as they 
appear at the top level of the spectral hierarchy. On the other hand, the analysis of a given 
system can be made as detailed as necessary (by using higher order approximations in each 
of the expansions (fT2j) b (JT3J), so the individual features are not overlooked behind broad 
universality. 

Work supported in part by the Sloan and Swartz Foundations. 
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